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Abstract. The charging of conducting or alternatively insulating dust grains in a 
supersonic plasma flow with a directed photon flux is studied by the particle-in-cell 
method. The electron emission modifies the surface charge distribution on a grain and 
in the surrounding plasma. The charge and potential distributions on and around a 
dust grain are studied for different photon fluxes and different angles of the incident 
unidirectional photons with respect to the plasma flow velocity vector. Continuous 
and pulsed radiations are considered. We show that photoemission allows controlling 
the charge on conducting grains, and discuss the charging of stationary and spinning 
insulating grains. Interactions between positively charged grains can be stronger 
than for negatively charged grains. The simulations are carried out in two spatial 
dimensions, treating ions and electrons as individual particles. 
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1. Introduction 

The charging of an object in a plasma is one of the basic problems in plasma physics. The 
understanding of this process is important in studies of interactions between the plasma 
and the object, or between many objects in a plasma. The question is particularly 
important in studies of dusty plasmas, where a number of charged dust grains can form 
ordered structures, such as dust clusters, strings, and crystals [1, 2, 3]. In dusty plasma 
experiments, grains are usually levitated in the sheath region above the electrode, and 
they are charged negatively due to the high mobility of electrons. Such grains can be 
exposed to an ion flow. A wake and characteristic regions of enhanced ion density (ion 
focus) are observed behind dust grains immersed in flowing plasmas [4, 5, 6, 7]. The 
ion focusing and the corresponding potential enhancement are more conspicuous for 
supersonic flows [7], and can lead to the alignment of grains in a direction of the flow 
[8, 9, 10, 11, 12]. Analogous problems can be formulated also for larger objects moving 
with respect to a plasma, such as spacecrafts or meteoroids [13, 14]. 

In a space environment, dust is often exposed to electromagnetic radiation [15]. 
Radiation can be directional or isotropic, either due to background radiation or 
scattering of directed light [16]. The situation is relevant not only for dust in space, 
but also for dusty surfaces of larger lunar bodies. In the latter case, the dust on the 
surface is charged by the plasma and directed solar radiation. It has been argued that 
the shadowing of light can lead to strong electric fields and transport the dust above 
the lunar surface [17]. In laboratory plasmas, the radiation is either due to the plasma 
glow or an external light source, thus it is either isotropic or directed, similarly to the 
space environment [18]. If the energy of incoming photons is larger than the work 
function of the dust surface material, the photoelectron current contributes to the net 
current to the grain and should be included in the charging analysis [3, 19, 20, 21]. In 
several respects the physics of this process resembles that for electron emissive probes 
[22]. The differences between the two physical processes are found in the mechanisms 
for the electron emissions and to some extent also in the velocity distributions of the 
emitted electrons. Photoemission will change the total charge on the dust grain and the 
surface charge distributions, and it can lead to new types of interactions between dust 
grains. Structures comprising positively charged dust grains in a plasma in the presence 
of UV radiation have been discussed theoretically [23, 24], and observed in experiments 
[25, 26]. 

A theory describing the charging of dust grains with photoemission in a self- 
consistent way is difficult to develop. In particular, photoelectrons can modify the 
plasma in the vicinity of dust grains. Several theoretical studies consider simplified 
models [24, 27, 28]. To study more realistic problems, one should employ numerical 
simulations, which can account for non-linear and other possible phenomena, and model 
the charging of a dust in plasma in a self-consistent way. By numerical simulations, it was 
demonstrated that the charge of the dust cloud in a plasma discharge can be modified by 
UV radiation [29]. The potential structures around a positively biased spacecraft were 
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also studied numerically [30] . In these works the electrons were treated as a Boltzmann 
distributed background. Neither of these studies considered the self-consistent charging 
of isolated objects. 

One of our intentions here is to provide a more realistic model by including electrons 
(photoelectrons in particular) in the analysis. In a recent communication, we have 
demonstrated that UV radiation allows for an accurate control of the charge on an 
isolated conducting grain [31]. We have also showed that photoelectrons can modify 
and polarize the surrounding plasma, enabling stronger interactions between positively 
charged conducting dust grains. 

In the present paper, we study numerically the charging of conducting or 
alternatively insulating dust grains in a supersonic plasma flow with a directed photon 
flux. We analyze the charge, density and potential distributions for different fluxes and 
energies of photons, and for different angles between the incoming unidirectional photons 
and the plasma flow velocity vector. Continuous as well as pulsed radiation is considered. 
Using a particle-in-cell (PIC) method, we simulate the entire charging process in a 
collisionless plasma. The simulations are carried out in two spatial dimensions, treating 
ions and electrons as individual particles. We consider unidirectional radiation, which 
is relevant for a dust in space exposed to the solar radiation, as for example lunar dust, 
or a laboratory experiment with an external radiation source. 

2. Numerical code 

We have modified the numerical particle-in-cell (PIC) code used in our previous 
studies [7, 32, 33], by including a photon flux and the photoelectric effect [31]. We 
consider collisionless plasmas in a two-dimensional system in Cartesian coordinates. 
Both electrons and ions are treated as individual particles, with the electron to ion 
temperature ratio Tg/Tj = 100, and Tg = 0.18 eV. The ion to electron mass ratio is 
mi/rrLe = 120 in most of the simulations. As a control case we analyze also results for 
a conducting grain with mi/me = 36720 (to represent Neon). The plasma density is 
n = 10^° m~^, and the plasma flow velocity is = 1.5 Cs, with Cg denoting the speed of 
sound. Because of the large thermal velocity of electrons, the plasma flow is represented 
solely by the ion drift. 

A circular dust grain of radius of i? = 0.375 in units of the electron Debye length 
Xoe is placed inside a simulation box of size of 50 x 50 X^e- The grain is assumed to be 
massive and immobile, except for the simulations of the spinning insulator. Initially, the 
grain is charged only by the collection of electrons and ions. For a perfectly insulating 
grain, a plasma particle hitting the dust grain surface remains at this position at all later 
times and contributes to the surface charge distribution. To model a small conductor 
in this work, the charge is redistributed equally on the dust grain surface at each time 
step. Such an algorithm is simple to use and is also found in other numerical studies 
[34] , but it does not account for the electric dipole moment on the conducting dust grain 
as induced by the anisotropic potential distribution in flowing plasmas. The equally 
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distributed surface charge will not necessarily cancel electric fields inside the grain, and 
thus the algorithm is not adequate for grains larger than the Debye length or for grains 
of shapes different than spherical (or circular). This algorithm is different form the 
one used in our previous studies of the charge distribution on larger dust grains, which 
enforced constant potential within the dust grain [7, 32]. The computational expenses 
of that algorithm were lengthy simulations and strict constraints on shapes and sizes of 
simulated dust grains. 

A directed photon flux is switched on after approximately 40 ion plasma periods 
Tj. At this time, we can assume that the surface charge distribution on a grain has 
reached a stationary level. The code is run typically up to 50 ion plasma periods. Three 
different angles between the incoming photons and the direction of the ion drift are 
considered: a = {0°, 90°, 180°}. For conducting grains, the simulated photon flux is 
^hu £ (0.25,2.5) X 10^^ m~^s""^. This together with photon energies Ehu of 4.8, 5.5 and 
7.2 eV, gives a photon power density H G (1.9,28.8) Wm~^. The work function W of 
the conducting dust grains is taken to he W = 4.5 eV, which is close to work functions 
of many metallic materials [24]. For insulating grains, the photon energies Eh^, are 10.3, 
11.0 and 12.7 eV. This, together with the photon fluxes as for the case of conducting 
grains, gives a photon power density of if G (4.5, 50.8) Wm~^. The work function of 
the insulating grain is taken as = 10 eV, which implies that photoelectrons will have 
the same energies as for the conducting case. 

When a photon hits the surface of the dust grain, a photoelectron of energy 
E = Ehu — W is produced at distance / = svAt from the dust grain surface, where 
s is an uniform random number s G (0, 1], At is the computational time step and v is 
the photoelectron speed. Photoelectron velocity vectors are uniformly distributed over 
the hemicircle and directed away from the dust grain surface, that is in accordance with 
Lambert's law. 

To investigate the stability of the surface charge distribution on insulating grains, we 
simulate also spinning grains. Instantaneous rotation of angles (3 = {1°, 5°, 10°}, as well 
as continuous rotation with the angular velocities Q = {7r/180, 7r/36, 7r/18, 7r/2, tt, 27r} 
in units of rad/r^ (corresponding to the grain rotation by angles of 1, 5, 10, 90, 180, and 
360 degrees within Tj, respectively) are considered. The rotation starts at approximately 
one ion plasma period after the onset of radiation. As a control case we also rotate the 
grain throughout the whole simulation. 

3. Numerical results 

The present section is in two parts. First we consider the charging of a conducting or 
alternatively insulating grain in the presence of continuous radiation. This problem is 
followed by the results from the simulations with pulsed radiation. 
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Table 1. The total charge qt on a conducting dust grain for dilFerent photon energies 
Ehi, and different photon fluxes (f>/i^ for a = 0°, averaged over lOri. The relative 
charge fluctuations Aqt are also shown. The total charge qt is normalized with the 

1/3 

unitary two-dimensional charge qo = e ['^o(3r>)] i where e is an elementary charge, 
and 7T,o(3£)) is the plasma density in the corresponding three-dimensional system. The 
unit of go is [qa] = C/m. 



Eh.^A.S eV Eh, 



= 5.5 eV 




0.0 

0.25 

0.50 

1.25 

2.50 



755 
163 
19 
251 
795 



4 
19 
173 
18 

8 



-755 
-168 



273 
1330 



12 



4 
17 
258 
18 
7 



3.1. Continuous radiation 

The charge on a conducting dust grain exposed to a continuous photon flux becomes 
more positive with the onset of radiation and saturates within one ion plasma period. 
Some of the resuhs for a conducting grain have been presented before [31], but we 
include them also in the present work for completeness. The saturation charge on 
a conducting grain, which is summarized in Table 1, depends on the flux density and 
photon energy. For a sufficiently high photon flux, the grain becomes positively charged. 
For low fluxes, the saturation charge does not depend signiflcantly on the photon energy. 
For higher fluxes, high energy photoelectrons lead to a more positive dust grain. The 
relative fluctuations of the charge are largest for the grain with the smallest charge. 
The absolute and relative charge fluctuations are smallest for the case without photo- 
emission. The results for the total charge for E^y = 7.2 eV are very similar to the case 
of Ehu = 5.5 eV, and therefore they are not presented in Table 1. 

The floating potential on a positively charged grain for two highest photon fluxes is 
shown in Table 2 together with the corresponding results from analytical calculations. 
The analytical results for the floating potential in Table 2 are calculated for a balance 
of the photoemission ihu, ion ii, and electron ie currents to the grain: ie = ii + ihu- 
For consistency, we restrict our analysis to a two-dimensional case. The photoemission 
current can for this case be expressed by [1]: 



where e > 0, and it is assumed that the photoelectric yield and photoemission efficiency 
equal unity. In the present two-dimensional model with unidirectional photons, Ah^, = 
2R, and ^hu{2D) = c{^hu{i.D)/cY^^ ■, with the physical dimension of [^hv(2D)\ = 
Subscripts (2D) and (3D) stand for two-dimensional and three-dimensional cases, 
respectively. The ion current to a plane surface segment with area Ai due to singly 
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Table 2. The floating potential on a grain for different photon energies Ehi, and 
different photon fiuxes (f>/i^ for a = 0°. The results from the simulations 'i'fi^ sim as 
well as from analytical calculations caic are shown. 

£^,1^=4.8 eV = 5.5 eV 





^fl, Sim 


^'/;. calc 


^/i, Sim 


^/i, calc 


(IQi^m-^s-i) 


(V) 


(V) 


(V) 


(V) 


1.25 


0.17 


0.17 


0.28 


0.27 


2.50 


0.16 


0.36 


0.56 


0.66 



charged ions drifting at supersonic speed w^, can be approximated by 

/ e\l/ \ 

ii = Aino{2D)Vdeexp -— , (2) 



where we define the ion cross section for supersonic ion flow as Ai = 2R. The ion current 
is consistent with the current to a probe for retarding fields [35], but we replaced the 
ion thermal velocity by Vd and neglected the numerical constant by assuming that ion 
velocities are unidirectional and normal to the probe surface at the sheath edge. We note 
that the ion current to the positively charged grain is negligible due to small thermal 
velocity of ions, but nevertheless we include it in the calculations for completeness. 
Since the grain radius R is comparable to the electron Debye length, we use a general 
expression for the orbit-motion-limited (OML) current to the conducting cylinder, to 
calculate the electron current to the grain [35]: 



1 r 



erf 



-7 



rl/R^ 



R 

H — exp(— 7) I 1 — erf 



where 7 = —e'^/kTe, = 27rr, and r^ is the sheath radius, which in our calculations 
is set to = 3-R. We introduced the error function as erf(x) = 2/^/7r Jq exp{—y'^)dy. 

The density and potential distributions around a conducting dust grain depend on 
the fiux and energy of the photons. For photon fiuxes of = 0.25 x 10^^ m~^s~^, 
when the grain is negatively charged, we observe an ion focusing in the wake [7]. The ion 
density in the focusing region is rij ^ 1.2noi, where noi is the undisturbed ion density 
far from the grain. This result is smaller than for the corresponding case without 
photoemission where we had n.j ~ 2.2noj. The ion focusing is destroyed for positively 
charged grains. In this case, ions are slowed down and defiected in front of the grain. 
Consequently, a region of an enhanced ion density is formed in front of the grain, while 
downstream from the grain there is a distinct boundary between the wake and the 
undisturbed plasma, see Fig. 1. The shape of the enhanced ion density region depends 
on a: it is more pronounced and located closer to the dust grain surface for a = 0°, and 
further from the grain for a = 180°. For a = 90°, an asymmetry in the enhanced ion 
density is observed [31]. 
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Table 3. The width w and length d of the ion wake behind a positively charged 
conducting dust grain for different photon energies Ehi, and different photon fluxes 
^hiy for a = 0° . The unit of w and d is the electron Debye length Xoe- The ion wake 
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— 9 — 1 
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The wake in the ion density behind a conducting grain (a region where rij < O.Snoi), 
scales with the photon flux and photon energy, increasing for increasing fluxes and 
energies. The measured spatial extent of the wake is summarized in Table 3. The ion 
wake corresponds to the white region behind the grain in Fig. 1. 

The potential around the positively charged conducting dust grain is polarized 
for higher photon fluxes. In Fig. 2, the potential distribution around the conducting 
dust grain is shown for different angles of incidence of photons with the flux = 
2.5 X 10^^ m"^s"^. The potential is negative behind, and positive in front of the grain. 
The polarization of the plasma is most conspicuous for a = 180°. 

The results from simulations with a more realistic ion mass are in accordance with 
the results obtained with the reduced ion mass. For a realistic ion mass, the total charge 
qt on a grain without photoemission is qt = — 1883go- This result is more negative than 
for simulations with reduced ion mass. The ratio of the saturation charges for different 
ion masses is qt,i/qt,2 = 2.5. It is close to the ratio 

Qo,i ^ In (7i/27r + l) ^ 
go,2 ln(72/27r + l) 

where indices 1,2 refer to different ion to electron mass ratios for ion masses rui = 
36720me, and = 120me respectively, and Qo is a theoretical charge on a grain in a 
stationary plasma in a two-dimensional system, given by 

In (5), Kq and Ki are modifled Bessel functions, R is the radius of a grain, 
7 = mi/me, and ^'j/ is a floating potential of the grain, here given by 

In (6) it is assumed that cold ions are reaching the surface of the large conducting 
object at the Bohm speed. A more detailed discussion on Equations (5) and (6) is given 
elsewhere [32]. 
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Figure 1. The ion density around a conducting dust grain exposed to the photon 
flux = 2.5 X 10^^ m^^s^^ of energy E/^, = 4.8 eV averaged over nine ion plasma 
periods r^. a = 0° and the plasma flow is in the positive x direction. The white region 
corresponds to ion densities below 0.5noi. 



Without photoemission, ions are streaming out of the ion focus with a wider angle 
for reahstic ion masses as compared to the case with a reduced ion mass. This is due to 
different ion drift velocities for the two cases. In both cases the ion drift is = 1.5 Cg, 
with the speed of sound given by Cs = \J i^(Te + 5Tj/3)/mj, in the plasma far from the 
grain. In the vicinity of the grain, the plasma parameters are modified due to particle 
trapping and sheath formation. With photoemission, the saturation charge and the 
wake are similar for both ion masses. The length of the wake is the same, while the 
width for larger ion masses is larger by 5%. The charge saturates within one ion plasma 
period for both cases. The ion plasma period for the realistic ion mass is approximately 
20 times larger than for the reduced ion mass. With another simulations for grains with 
the radius of 2R, we find that the saturation charge is approximately twice the charge 
value for the grain with radius of R. 

The charging of an insulating grain exposed to a continuous radiation differs from 
the conducting case. The saturation in the charging characteristics is observed for 
photon fluxes ^hu = 0.25 x 10^^m~^s"^ and = 0.50 x 10^^m~^s~\ when the total 
charge on the dust grain remains negative. For photon fluxes and energies high enough 
to change the sign of the total charge on the grain, the charge does not saturate within 
the time-span of our simulations. In all cases, the charging depends on the angle of 
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Figure 2. The potential around a conducting dust grain exposed to the photon flux 
= 2.5 X 10^3 m-2s-i of energy Eh^ = 4.8 eV for a = 0° (a), a = 90° (b), and 
a = 180° (c). The data were averaged over a time interval of nine ion plasma periods 
Ti. The plasma flows in the positive x direction. 

incidence, see Fig. 3. For lower fluxes the charge is getting less negative with increasing 
a. For higher fluxes, the charge can become positive, and then negative again within a 
few ion plasma periods. This is not the case for a = 180°, for which the charge increases 
towards more positive values. 

With the onset of the photon flux, we observe the development of an electric dipole 
moment on the grain which is antiparallel to the direction of the incident photons, see 
Fig. 4. This electric dipole moment due to the photoelectrons does not saturate for high 
photon fluxes, and it is stronger than the electric dipole moment induced by the ion 
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Figure 3. The total charge on an insulating dust grain as a function of time for 
different photon fluxes and angles of photon incidence a. Squares correspond to the 
photon flux ^hu = 2.5 x 10^^ m^^g-i^ triangles to <^hu = 0.5 x 10^^ m^^g-i j,^^ 
photon energy is E}i^ = 11.0 cV. The results arc smoothed with a moving box average 
filter for presentation. 



The density and potential distributions around an insulating grain evolve in time, 
see Figs. 4 and 5. The potential distribution is influenced by the electric dipole moment 
due to photoemission. This moment becomes smaller when the total charge is negative. 
For a positively charged grain, the ion focusing region in the wake is destroyed, and the 
wakes behind the dust grain with a = 0° and a = 180° are similar to the conducting 
case. The wake is strongly asymmetric for a = 90°. For high photon fluxes, when the 
charge on the grain reaches negative values, the wake behind insulator becomes smaller, 
and the ion focus can be retrieved. The asymmetric charge distribution for a = 90° is 
present also after the closure of the wake, as shown in Fig. 5b). 

In Fig. 6, we illustrate the difference 6 between the density of Boltzmann distributed 
electrons that would correspond to the calculated potential and the actual electron 
density: 6 = rzeo exp [e^'/ fcTg] —ne, where e > is the magnitude of the electron charge. 
Results for both conducting and insulating stationary grains are shown in Fig. 6. Before 
the onset of the photon flux the electrons can be well approximated by the Boltzmann 
distribution. With photoemission, the electrons are no longer Boltzmann distributed. 
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Figure 4. The potential around an insulating dust grain exposed to the photon flux 
= 2.5 X 10^^ m-^s-i, a = 0° (a) and a = 90° (b) of energy Eh^ = 11.0 eV 
averaged over two ion plasma periods: t € (39.5, 41.5)Ti (top) and t E (48.0, 50.0)ri 
(bottom). The plasma flows in the positive x direction. 

The largest discrepancies for conductors are associated with a surplus of electrons due 
to the photoelectron emission, and to a region of an enhanced ion density in front of 
the dust grain, where electrons are underrepresented. For insulators the electric dipole 
governs the potential in vicinity of the grain. 

Instantaneous rotation of the insulating dust grain with an angle f3 has little effect 
on the grain charging characteristic. For a = 180°, no significant change is observed 
in the potential and density distributions, while for a = 0° the rotation leads to weak 
asymmetries there. The asymmetries are more pronounced for larger (3. For a = 0° the 
charging characteristics are similar to the case without rotation, but the charge becomes 
more negative at a slightly slower rate with increasing f]. 

Continuous rotation by an angle of ir within one ion plasma period significantly 
modifies the charging of the grain, see Fig. 7 (a). The rotation of a grain redistributes 
the charge on the dust grain surface, and lowers the total charge on the grain. After 
arresting the grain rotation, the charge becomes more negative for a = 0° and a = 90°, 
while for a = 180° it can saturate with a quadrupole moment in the surface charge 
distribution. 
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Figure 5. The ion density around an insulating dust grain exposed to the photon 
flux = 2.5 X lO^^ m-2s-\ a = 0° (a) and a = 90° (b) of energy E,,,, = 11.0 eV 
averaged over two ion plasma periods: t g (39.5, 41.5)ri (top) and t e (48.0, 50.0)ri 
(bottom). The plasma flows in the positive x direction. White regions correspond to 
ion density levels below 0.5noi. 

With steady state rotation, the total charge oscillates in time with the mean charge 
value lower than on the conducting grain with the corresponding parameters for the 
photon flux. The period of oscillations depends on the angular velocity of the grain, see 
Fig. 7b). The oscillations are not observed for very slow angular velocities. There is 
little difference in the grain charging characteristics for different starts of the rotation 
of the grain. For a grain spinning throughout the whole simulation, the total charge 
before the onset of photoemission is less negative than on a stationary grain. This is 
due to the charge redistribution, which prevents the development of a strong electric 
dipole moment. However, with photoionization, the charging characteristics are similar 
to the case when the dust grain starts spinning after the radiation onset. The charge 
redistribution on spinning grains tilts the electric dipole moment on insulating grains. 
The strength of the electric dipole moment oscillates together with the total charge on 
the grain. Simultaneously, the wake becomes asymmetric and its size oscillates in time. 
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Figure 6. The difFcrcncc S between the density of Boltzmann distributed electrons 
that would correspond to the calculated potential and the actual electron density is 
shown for the case with (solid line) and without (dashed line) photoemission. Both 
conducting (a) and insulating (b) dust grains are considered for a = 0°. 



3.2. Pulsed radiation 

The charge on the conducting grain exposed to a radiation pulse is more positive during 
the illumination. After the pulse, the charge recovers to the previous value (before the 
pulse) within approximately one ion plasma period. The charge recovery is initially 
fast and then continues at a slower rate. Initially, the charge recovery can be well 
approximated by an exponential function of the form q = goexp[— t/r], with the time 
constant r = 3.45 x 10~^s for = 2.5 x 10^^ m~^s~^, and r = 4.56 x 10~^s for 
^hu = 0.5 X 10^^ m~^s~^. These time constants are comparable with the electron 
plasma period Tg = 3.53 x 10~^s, which suggests that initially the charge recovery is 
primarily due to electrons. The time constant for = 0.5 x 10^^ m~^s~^ is larger than 
Te because the maximum charge is close to zero in this case, and the ions contribute 
initially to the charge recovery. After a time interval of 2re, the time constant r is 
larger and reaches r ^ 0.5rj at the end of the recovery for both cases. The charging 
is shown in Fig. 8a) together with points corresponding to the exponential fits. For 
clarity of presentation, we do not show continuous exponential fits, but only regularly 
spaced points corresponding to the locally fitted curves. Approximately one ion plasma 
period after the switch off, a small overshoot in the charging characteristic is observed 
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Figure 7. The total charge on an insulating dust grain rotating with an angle of tt over 
one ion plasma period (a) and continuously rotating after the start of photoemission 
(b). Continuous radiation with <f>;i^ = 1.25 x 10^^ ni~^s~-'^ (squares correspond to 
^hv = 2.5 X 10^^ ni~^s~"'^) is switched on at i = SOr^. Triangles correspond to the 
dust grain spinning throughout the whole simulation, a = 0° for = tt, a = 90° for 
n ~ 0.57r, and a = 180° for = 27r in units of t^^ . The results are smoothed with a 
moving box average filter for presentation. 



for higher photon energies. The charging at given photon fluxes depends only little on 
the photon incident angles a. 

The charging after a series of three pulses is similar to what is found for a single 
pulse with the relevant photon flux and energy. Each pulse corresponds to a peak in 
the charging characteristics in Fig. 8b). The height of each peak does not change much 
with the time interval between the pulses. The trough is less negative for time intervals 
between the pulses that are shorter than the charge recovery time. 

The electrostatic potential around the conducting dust grain exposed to radiation 
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Figure 8. The total charge on a conducting dust grain exposed to a single radiation 
pulse (a) and to a pulse series with different time intervals between pulses Atp (b). 
In plot (a) different symbols represent regularly spaced data points corresponding to 
local exponential fits with different dime constants r. In plot (b) triangles correspond 
to = 0.5 X 10"^^ m~^s~^ and squares to = 2.5 x 10^^ m~^s~^. In both cases 
a = 0°, and Eji^, ~ 5.5 cV. The results are smoothed with a moving box average filter 
for presentation. 



pulses is polarized as in the case of continuous radiation. During the pulses, the potential 
behind the dust grain is negative, and resembles the case of the conducting grain with 
continuous radiation, see Fig. 2. This region remains negatively charged also between 
the pulses. Within a time interval of approximately 1.5rj after the last pulse, the positive 
potential region in the grain wake is rebuilt: first in the vicinity, and then further away 
from the grain. At the same time, the region with net negative charge becomes less 
pronounced and moves further downstream from the dust grain, slower than the ion 
drift speed. During the pulses, the wake potential in the vicinity of the grain oscillates 
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Figure 9. The total charge on an insulating dust grain exposed to the pulsed radiation 
as a function of time for different time intervals between pulses Atp and different photon 
fluxes: = 0.5 x 10^^ m'^s"^ (a) and ^hi, = 2.5 x 10^^ m^^s"^ (b). Squares 

correspond to a = 180°, triangles to a = 0°. The photon energy is E^i, ~ 11.0 eV. 
The results are smoothed with a moving box average filter for presentation. 



with the frequency of the pulses, see Fig. 11. These oscillations propagate into the wake, 
but are heavily damped further away from the dust grain, and diminish after the last 
pulse. 

The ion density behind a conducting dust grain is being enhanced between the 
pulses. This process is terminated by the start of a successive pulse, and the enhanced 
ion density regions move away from the dust grain in the ion drift direction, but at 
a lower speed. The full recovery of the ion focus occurs after approximately one ion 
plasma period from the last pulse. The ion density wake in front of the grain is located 
closer to the grain for a = 180° than for a = 0°. There is also a depletion in the electron 
density in the region corresponding to the ion wake originating from a positively charged 




X [units of Xoe] 



Figure 10. The potential around an insulating dust grain after a series of pulses 
with Atp = l.Or,. = 2.5 x 10^^ m^^s-i, E,,^ = 11.0 eV, and a = 180°. The 

minimum in the potential is ^' = —12.2 in units of kTe/e on the surface of the dust 
grain. Potentials lower than ^' = — 5 kTp/e are coloured black. 

grain. After the switch-off, photoelectrons are rapidly redistributed, but the depletion 
in the region corresponding to the wake remains until the wake is filled with ions. 

An insulating dust grain exposed to a single radiation pulse has charging 
characteristics similar to the conductor only for low photon fluxes and a = 0°, when the 
total grain charge remains negative, see Fig. 9a). For the low photon flux and a = 180°, 
the total charge on a dust grain is low during the pulses. In this case, the charge does 
not recover within one ion plasma period after the pulse, but it reaches half of the charge 
value before the onset of radiation. Longer simulations show that the recovery time is 
approximately 20 ion plasma periods in this case. 

For high photon fluxes, the charge during the pulse is positive, see Fig. 9b). After 
a single pulse with a = 0° the charge is more negative than before the pulse, and less 
negative when a = 180°. The total charge in the trough between subsequent pulses is 
getting more negative for a = 0°, while for a = 180° it is more negative for longer time 
intervals between pulses, and less negative for shorter intervals. In all insulating cases, 
the total charge after a series of pulses can be more negative than after a corresponding 
single pulse, and the full charge recovery takes usually several plasma periods. 

The radiation pulses modify the potential and density patterns in the vicinity of 
the grain. For photons with a = 0°, we observe an enhancement in the electric dipole 
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Figure 11. The potential variations at different distances Ax from the rear of the 
conducting grain for y ~ 25.8 in units of Xjj^ exposed to the radiation pulses as a 
function of time. = 2.5 x lO^^ m-^g-i, Eh^ = 5.0 eV, and a = 180°. 

moment on the dust grain. For photons with a = 180°, the positive surface charge 
is accumulated on the front and rear sides of the dust grain. During such pulses, the 
electric dipole moment is parallel to the ion flow (antiparallel to the photons direction), 
and after the pulses, the quadrupole moment in the surface charge distribution develops, 
see Fig. 10. The quadrupole moment diminishes in time, faster for low photon fluxes, 
and the dipole moment in the surface charge distribution as well as the electrostatic 
potential distribution are recovered. 

In the wake, the ion focusing is rebuilt behind the dust grain after the radiation 
pulses. For a = 0°, this region is similar to the one before the pulses, while for a = 180°, 
at the same time instances it is wider and weaker for low fluxes, and spatially narrower 
with stronger focusing for high fluxes. Electrons are Boltzmann distributed after the 
pulses. 

4. Discussion 

Photoemission provides an electron source on the irradiated side of the grain and 
modifies the dust grain charge. For sufficiently high photon fluxes, the charge on 
the conducting grain becomes positive and saturates within one ion plasma period. 
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Positively charged conducting grains slow down and deflect flowing ions. As a result, 
a region of enhanced ion density forms in front of the grain, while behind the grain 
a substantial wake in the ion density is formed, see again Fig. 1. The wake in the 
ion density scales with the photon energy and flux, being larger for higher fluxes and 
energies. Hence, the wake size is proportional to the charge on the grain. Photoelectrons 
with energies higher or comparable to the electron thermal velocity can easily be lost on 
the grain surface, while with higher energies they are more likely to escape the trapping 
potential of the grain. This together with the photoemission rate, which is proportional 
to the flux, explains the development of more positive charge on the grain for high 
energetic photons and high fluxes [31]. The angle between incoming photons and plasma 
flow direction has little effect on the potential distribution around conducting grains. 
Photoelectrons contribute in neutralizing enhanced ion density regions. For this reason, 
the region of enhanced ion density is located closer to the front of the grain for a = 0°, 
when the grain charge can be effectively shielded by the photoelectrons, and further away 
for a = 180°, when the photoelectrons are produced on the shadow side. The electrons 
penetrate into the ion wake, due to their high mobility. The resulting imbalance between 
ion and electron densities in the vicinity of the grain leads to polarization of the plasma, 
see again Fig. 2. This allows for strong interactions between many positively charged 
grains in flowing plasmas. The electrons are no longer Boltzmann distributed. 

To calculate the charge on positively charged conducting grains, we use a two- 
dimensional capacitance model [32]. In (5) we use the simulation results for the floating 
potential For photon fluxes of = 1-25 x 10^^ m~^s~^, the calculated charge is 
qt ^ 250go for both photon energies. This result is close to the data shown in Table 
1. For ^hu = 2.5 X 10^^ m"^s~^ we have qt ~ 410go for Ehu = 4.8 eV, and qt ~ 732go 
for Ehi, = 5.5 eV, which is lower than the simulation results. However, if in equation 
(5) we formally substitute by X^e for $^2^ = 2.5 x 10^^ m~^s~^, then the results 
are close to the simulation results: qt ~ 734go for Ehu = 4.8 eV, and qt ~ 1312go for 
Ehu = 5.5 eV. This suggests that for lower photon fluxes, and low positive potentials 
of the grain, the ions can effectively contribute to the shielding of the grain, while for 
more positive potentials, the grain potential is predominantly shielded by electrons. 

The analytical solution for the floating potential is in a good agreement with the 
simulation results for low energy photons. For high energy photons the analytical 
calculations give more positive potentials than obtained from the simulations. This 
is due to thermalization of the photoelectrons. The temperature of low energy photo- 
electrons is higher than the mean temperature of the background electrons, but the 
corresponding velocities are still within the thermal spread of the background electrons. 
High energy photoelectrons are effectively slowed down by the grain and interact with 
background electrons. We flnd that the average rate of deceleration of the high energy 
photoelectrons is 25% of their initial energy. With this correction for high energy 
photons, the analytical calculations for the floating potential give values close to the 
numerical results. 

An electric dipole moment develops on insulating grains due to the photoemission. 
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It is oriented antiparallel to the photon direction. Neither the electric dipole moment 
nor the charge saturate on such grains within the simulation time. For a = 180° 
the charge is more positive in time, while for other angles of incidence it recovers to 
negative values. For a = 180° both rear and front of the grain are positively charged. 
The depletion in the ion density behind the grain does not allow electrons to neutralize 
the rear charge. With increasing charge on the rear of the grain, a weak quadrupole 
moment develops (with negatively charged grain sides tangential to the flow), and the 
total charge increases towards positive values. For other angles, the non-irradiated side 
of the grain becomes more negative when the photoelectric current exceeds the electron 
current. Therefore, a depletion of the wake and a recovery of the ion focusing region 
are observed, see again Figs. 4 and 5. The surface charge distribution on insulating 
grains for a = 90° leads to more pronounced asymmetries in the ion density than on 
conducting grains. The electrons are no longer Boltzmann distributed also in the case of 
insulators, but a dominant contribution here is due to the electric dipole development. 
Photoelectrons are absorbed by the high positive charge on the dust grain, and they 
neutralize enhanced ion density regions in front of the grain. 

Rotation of the insulating dust grain redistributes the charge on the grain surface. 
Without photoemission but retaining the directed ion flow, the total charge on fast 
spinning grains becomes less negative. The electric dipole moment on the surface 
diminishes and the value of the total charge is similar to the conducting case. With 
photoemission, the electric dipole moment on the spinning dust grain is still present for 
the photon fluxes considered in this work. It is tilted by an angle with respect to the 
direction of radiation, and the positive charge on the irradiated side of the dust grain is 
being neutralized when it reaches the shadow side. Due to the depletion of the wake in 
the ion density, the total charge oscillates on fast spinning grains, see Fig. 11. With the 
spinning grain, the symmetry in the ion wake is destroyed near the grain surface. For 
sufficiently fast rotation, the redistribution of the negative charge bends ion trajectories 
and leads to wake erosion. This process continues until the ion density is rebuilt in the 
vicinity of the grain and the region of reduced density is detached from the grain. The 
electron density in the wake increases, and so does the total electron current to the grain. 
When the charge becomes less positive, and the electron current to the grain decreases, 
photoemission leads to the formation of the new wake in the ion density. The closing of 
the wake and the resulting oscillating total charge will occur only if the erosion of the 
ion wake is substantial. If the rotation of the grain is slow and the photoemission rate 
is high, the wake will not close and detach. 

The present discussion considers grains with spherical, cylindrical and oblate 
shapes. Irregularities on the dust grain surface can lead to variations in the surface 
charge distribution and the wake can be perturbed also for smaller angular velocities of 
the grain. It is also noted, that spherical grains with inhomogeneous surface properties 
will spin due to angular momentum transfer from ions [36]. Such inhomogeneities can 
also lead to complicated surface charge distributions, as illustrated elsewhere [32]. 

Due to the high inertia of the dust grains, the rotation of the grains will be slow 
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in most experiments [37]. Hence, irregularities in the dust grain surface that give rise 
to a comphcated surface charge distribution will be the main factor for the charge 
saturation on the insulating grain in the presence of directed radiation. This will be 
valid also for slowly spinning grains. On the other hand, the surface charge distribution 
on the insulating dust grain exposed to directed radiation will lead to strong electric 
fields within the grain. This effect will be more pronounced on grains with surface 
irregularities. The irregularities will eventually be destroyed by strong electric fields. 
This will be valid also for stationary plasmas and is similar to the sterilization and 
destruction of bacteria by means of plasma used as a source of UV radiation [38, 39]. 

Photoemission provides a method for controlling the charge on conducting grains 
both in vacuum [18] and in plasma [31]. For perfectly insulating grains, the total charge 
saturates only when the total grain charge remains negative, but it depends on the 
photon incidence angle. The charge should also saturate for slow spinning insulators 
due to surface irregularities. It oscillates for fast spinning insulating grains, however. 

After pulses of radiation, the charge recovers to the value from before the onset 
of radiation. The recovery is initially mainly due to electrons, and then both electrons 
and ions. The charging can be approximated by an exponential function with the time 
constant r that is initially comparable with the electron plasma period and then larger. 
This is in agreement with the previous results from dust grain charging simulations [32]. 
The charge on conducting grains recovers within one ion plasma period. The small 
overshoot in the charging characteristics at approximately one ion plasma period after 
the switch off can be attributed to the ion response due to the reduced ion mass [32], 
and the formation of the ion focus. For insulators, the charge recovery may take up to 
20 ion plasma periods. This is due to a complicated surface charge distribution on the 
grain. For a = 180°, the positive charge on the rear of the grain is reduced slowly after 
the pulse because of the reduced density in the wake, while the front side of the grain 
is positively charged by the ion flow. A quadrupole moment is present in the surface 
charge distribution for several ion plasma periods after a pulse. Consistently, the charge 
between and after the pulses is less negative for a = 180° than for other angles. For 
other angles, the surface charge distribution leads to a more negative total charge after 
the switch off as compared to the dust grain charge previous to photoemission. 

The photon fluxes considered for conductors in this work, can be achieved by 
commercially produced sources of UV radiation (e.g., low pressure mercury lamps) [38]. 
In case of lamps it would be necessary to coUimate the light, but in case of UV lasers 
the energies of E^i, G (4.8, 7.2)eV, corresponding to A G (172,258)nm, can be achieved 
for instance by excimer lasers used in photolithography [40]. The photon energies in 
the range E^,^ G (10.3, 12.7)eV, corresponding to A G (97, 120)nm, are more difficult to 
obtain, but can still be achieved for instance by free electron lasers [41]. In our work 
we assumed the work function for the insulator to he W = 10 eV, which is similar to 
values from experiments with ice [42, 43]. However, insulating grains can have lower 
work functions. The work function of pure ice is W ^ 8.7 eV, and can be significantly 
lower if the ice contains impurities, as it is expected in the atmosphere [44]. Sodium 
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silicate glass can have work function as low a.s W = 6 eV [45]. Therefore, on such 
insulators similar effects, as demonstrated in this study, should be achieved by lower 
photon energies. 

By illuminating a grain using short pulsed lasers, it is possible to modify the charge 
on the grain and excite potential oscillations in the wake. The other dust grain, if 
located in the wake, will experience oscillations. Since the charge on the illuminated 
grain is determined by photoemission, the motion of the particle would provide non- 
invasive diagnostics for measuring the charge on the grain located in the wake of the 
illuminated one. On the other hand, continuous illumination of the grain placed in the 
wake will fix the grain charge and allow for accurate study of the wake of other grain. 

The results presented here does not depend on the ion to electron mass ratio, except 
for the saturation charge on the dust grain without photoemission and the charging rate. 
However, the charge on the grain with photoemission does not change with the ion to 
electron mass ratio. This is because the photoemission current does not change with the 
ion to electron mass ratio, and the ion current to positively charged grains is negligible. 

There are certain limitations of our model. The radiation is assumed to be 
unidirectional and photoelectrons to be monoenergetic. In many problems in laboratory 
plasmas, UV radiation will be more isotropic due to the plasma glow and light scattering, 
while photoelectron energies will be statistically distributed. Isotropic radiation will 
cause more homogeneous distributions of the photoelectons and the grain surface 
charge. We considered perfectly insulating and perfectly conducting dust grains. Finite 
conductivity due to impurities and resistivity can modify the results, especially for 
insulators. These issues were not considered in this work. 

5. Conclusions 

The results from numerical simulations of charging of isolated dust grains in flowing 
plasmas with photoemission were presented for perfectly insulating and perfectly 
conducting grains. By means of photoemission, the total charge on a conducting 
grain can be effectively controlled. The charge control on insulating grains is more 
difficult since no charge saturation is observed for stationary grains, and the charging 
characteristics depend on the angle between the incident photons and the plasma flow. 
For insulating grains, surface charge irregularities and rotation of the dust grain can 
redistribute the surface charge on the grain. Fast spinning of the grain results in 
oscillations of the value of the total grain charge and the density wake behind it. 

During photoemission, the electrons are non-Boltzmann distributed in the vicinity 
of the grain. This makes a theoretical analysis of the problem difficult. The plasma is 
polarized in the vicinity of the grain, which can give rise to strong interactions between 
many dust grains. For insulators the interactions are controlled by a strong electric 
dipole moment on the surface, antiparallel to the direction of radiation. After pulses of 
radiation, the charge, density and potential distributions recover to the conditions from 
before the photoemission. The recovery takes approximately one ion plasma period for 
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conducting grains, and several times longer for insulating grains, due to the complicated 
charge distributions on the dust grain surface. 

By a fine adjustment of the charge with the use of photoemission, the coagulation of 
the dust grains can be induced due to large relative fluctuations of the charge when the 
total charge on a grain is small. Both continuous and pulsed radiation should allow for 
non-invasive diagnostics of the charge and wake structure in dusty plasma experiments. 
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